Half-quantum vortex and d-soliton in Sr 2 Ru04 
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Assuming that the superconductivity in Sr^RuCU is described by a planar p-wave order parame- 
^■y ter, we consider possible topological defects in S^RuCU. In particular, it is shown that both of the 

d-soliton and half-quantum vortex can be created in the presence of the magnetic field parallel to 
the a-b plane. We discuss how one can detect the d-soliton and half-quantum vortex experimentally. 

G> ' PACS numbers: 74.20.-z, 74.80.-g, 74.25.-q 

(N : 

It has been suggested that the unconventional superconducting state of Sr2Ru04 is described by the planar spin- 
triplet p-wave order parameter with broken time reversal symmetry in analogy to the 3 He A-phase. As known from 
the example of superfluid 3 He, one of the hallmarks of the triplet superconductivity is the presence of a manifold 
of topological defects. Thus, we expect that the creation and detection of topological defects in Sr2Ru04 (or 
spin-triplet superconductor) will provide further insights about the nature of the unconventional superconducting 
I 1 state of Sr2RuC>4. In this sense, the study of the topological defects in the planar p-wave superconducting state with 
^ | broken time reversal symmetry is of great interest. 

More specifically, it has been proposed [HJJ] tna t the superconducting order parameter in this system is described 
| by the planar p-wave form written as, 

S A af3 {k) = d{k) ■ (aio- 2 ) a p (1) 

, with 

§; d(k) = Ad(k ± ik 2 ) , (2) 

where cr M (p = 1,2,3) are Pauli matrices and a, (3 represent spin f or J,. Here kj(j = 1,2) represent the projection of 
the unit wave vector k along two perpendicular directions e\ and e.2 in two dimensional space. This order parameter 
describes the Cooper pair state with the zero spin projection on d and with the unique projection of the pair orbital 
angular momentum given by I = e\ x e.2- In S^RuC^, due to the spin-orbit coupling, d is forced to be parallel to ±c 
and k is the quasi-particle momentum in the a-b plane. 

Indeed the spontaneous magnetization seen by muon spin relaxation experiment [^| and flat Knight shift seen by 
NMR § seem to be consistent with this picture. On the other hand, the origin of the spontaneous magnetization 
seen by muon spin relaxation experiment is somewhat mysterious since we do not expect such a magnetization in a 
^ • homogeneous system. 

It is important to notice that the superconducting ground state described by the order parameter of Eq.|| is doubly 
degenerate. We can designate these two ground states by the angular momentum l z = ±1, where 1 is parallel to the 
c-axis. Sigrist and Agterberg Q proposed recently that there will be in general a domain wall between l z = 1 and 
l z = — 1 states, which we shall call [-soliton |s) in analogy to the case of superfluid 3 He-A. It is likely that such a 
soliton is magnetically active, so it may be an origin of the spontaneous magnetization seen in muon spin relaxation 
experiment || . In particular, in a magnetic field H \ | c, only one of these degenerate states is favored Jj|[l0| . Therefore, 
it is possible to control f-solitons by a magnetic field parallel to the c-axis. jjj They also proposed that these f-solitons 
would provide very efficient barriers for the vortex motion and this effect is possibly related to the pinning of vortices 
observed in Sr2Ru04 below T — 30mK pl| . However, in this experiment the magnetic field is applied in a direction 
perpendicular to the c-axis. As we will see later, so-called d-solitons appear to be more appropriate than Z-solitons in 
this configuration. Also, for the Z-solitons, it is rather difficult to estimate the soliton energy and to make a further 
quantitative prediction. 

The purpose of this paper is to propose an alternative model for the appearance of the spontaneous magnetization 
and the mechanism of the pinning of vortices; <i-soliton and half-quantum vortex. The d-soliton is a domain wall 
between d \\ c and d || — c as in superfluid 3 He-A. We believe that d is parallel and antiparallcl to c, because they are 
forced to be parallel to the angular momentum / (or —I) due to the spin-orbit coupling characterized by an energy 
scale fid fl§]- Therefore, if we use Qd as a parameter, we can calculate the energy and shape of the e?-soliton provided 
that fid <C A(T), where A(T) is the superconducting gap. Unfortunately we do not know the precise value of f^, 
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but it may be about ^jA(T). In this picture, moving d-soliton generates the local magnetization which can result in 
the spontaneous magnetization seen in muon spin relaxation experiment. One can also generate a large number of 
d-solitons by applying a burst of high frequency microwave with frequency ~ f2 d sent parallel to the a-b plane. J8| 

As in superfluid 3 He-A, each d-soliton is terminated by a pair of half-quantum vortices. Jl3|] We find that these 
pairs of half-quantum vortices are more stable than the usual single quantum vortex in the superconducting state in 
the presence of the magnetic field parallel to the a-b plane. This means that the usual single quantum vortex would 
spilt into a pair of half-quantum vortices connected by the d-soliton. PlJlifl In this case, these objects would provide 
an extremely efficient pinning mechanism of vortices in Sr 2 Ru04. Also the half-quantum vortices should exhibit a 
clear electron spin resonance (ESR) signature. Further we believe that these objects are visible by the scanning 
tunneling microscopy (STM) imaging and by micromagnetometry developed by Kirtley et al [ fL5| used in high T c 
cuprate compounds. 

Free energy of the conventional single vortex when the magnetic field is parallel to the a-b plane 

Let us assume that the magnetic field is parallel to the a-axis. Then the free energy of the conventional vortex with 
the flux quantum fa — hc/2e is obtained within the London approximations as 

where £ is the coherence length and A is the magnetic penetration depth. The magnetic penetration depth A is 
related to the superfluid density p s {T) by A~ 2 = ^rp a (T). When the magnetic field is parallel to the a-axis, in the 
anisotropic system like Sr 2 Ru04, A and £ should be reinterpreted as A = VA&A C and £ = V£b£c- Here Af, jC and £b,c 
are the magnetic penetration depth and coherence length in the b and c directions respectively. 

d-soliton and a pair of half-quantum vortices 

There exists huge anisotropy in the in-plane and out-of-plane transport properties in Sr2Ru04. Thus Sr2Ru04 may 
be regarded as an effectively two dimensional system. The large anisotropy or the effective two-dimensionality of the 
system forces the angular momentum of the Cooper pair to be parallel or antiparallel to the c-axis. In the p-wave 
superconducting state described by the order parameter described by Eq.jl] and EqJ|, the d vector is oriented along 
±Z in the presence of the spin-orbit coupling. Here we consider the case that the angular momentum I is uniform 
in the entire system. We can assume, without loss of generality, that I || c. We are interested in the deformation of 
the d configuration from the uniform case; for example, d || c. Any deviation from the uniform state would cost the 
energy associated with the spin-orbit coupling characterized by an energy scale fid fl2fl . However, we will show that 
the so-called d-soliton (a particular form of the d configuration) with a pair of half-quantum vortices can have lower 
energy than the conventional single vortex. Thus it is easier to excite a pair of half-quantum vortices with a d-soliton 
compared to single conventional vortex. In particular, a magnetic field parallel to the a-b plane generates very likely 
pairs of half-quantum vortices rather than usual vortices when the formers are stable. 

We consider the d-soliton that is a topological planar defect in the d configuration. The orientation of d changes by 
7r across the planar defect while d vectors at far distances are still along the c-axis. Typical configurations of d-soliton 
in the y-z plane can be found in Fig.l and Fig. 2 which we will explain later. We take y and z as the coordinates along 
6-axis and c-axis respectively. Now let us attach a pair of half-quantum vortices to the end points of the d-soliton of 
length R in y-z plane. In the case of an isolated half-quantum vortex, we have e i7r = — 1 factor in the order parameter 
due to phase winding around the half-quantum vortex. Therefore, an isolated half-quantum vortex cannot occur. On 
the other hand, if the half-quantum vortex is attached to the end points of the d-soliton, the disgyration in d at the 
same point compensates the phase tt so that there is no net change in the overall phase of the order parameter. 

In order to show that a pair of half-quantum vortices with the d-soliton is a lower energy excitation compared to 
single conventional vortex, we have to compare the free energies of two cases. The free energy required to create the 
d-soliton is obtained from 

h = \xnC 2 f d 3 r[£ \didtf + C 2 (l ~ dl)] , (4) 

ij 

where xn is the spin susceptibility, £d(T) = C(T)/Q.d(T) where C(T) is the spin wave velocity, and fid (T) is the 
longitudinal spin resonance frequency. Jl2fl 

On the other hand, d vector of the d-soliton can be parametrized by the following expression. 

d = cos ipz + sin ipy , (5) 
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where 



,/ x 1 ( z + R/2 z-R/2\ 

ip(y, z) = - arctan arctan , (6) 

2 V V V J 

where we put two half-quantum vortices at (y, z) — (0, R/2) and (0, —R/2). 

In the past, similar form of ip was also discussed in a different context in reg ard to 3 He. || As one can see, there 
is a discontinuity in ip across the line defined by —R/2 < z < R/2 and y = 0. The spatial configuration of the 
corresponding d around a pair of half-quantum vortices is shown in Fig.l. As one can see from the figure, the planar 
defect is parallel to the z-direction or c-axis. One can also consider the planar defect lying along the y-axis given by 

If y + R/2 y-R/2 \ 
ip = — I arctan arctan I , (7) 

where two half-quantum vortices are located at (y, z) = (R/2,0) and (— R/2,0). The configuration of the d vector 
using the above ip is shown in Fig. 2. One can easily see that the free energies, fd, associated with two possible d 
configurations are the same. 

The total free energy of the d-soliton and a pair of half-quantum vortices is given by 

fpair = \xnC 2 j dydz [K(V<S>) 2 i + £ \dj 3 \ 2 + £f sin 2 V] 

ij 

= l XN C 2 (7rK\n± +I 1+ I 2 ) , (8) 

where $ represents the phase of the order parameter which couples to the external electromagnetic field. The 
parameter K is defined by 

K _ Ps = 1 + 1/3*1 1 + l/3Ff (1 - g°) 

Psp l + l/3Ff l + l/3Fi(l-pg) ' 1 ' 

where p s and p sp are the superfluid density and the spin superfluid density respectively. F± and are the Landau 
Parameters and p® (= 1 — Y(T) and Y(T) is the Yosida function) is the superfluid density without the Fermi liquid 
correction. Q Notice that K{T C ) — 1 at T — T c and K{Q) — \^_\^p\ at T = 0. The temperature dependence of the 
parameter K is shown in Fig. 3 assuming that Fi = 9 and F° = 0. This choice of the parameters will be explained 
later. 

The first term in the first and second lines of Eq.|^ is the contribution from two half-quantum vortices and represents 
the fact that these half-quantum vortices repel each other. ii,2 are the contributions from the second and third terms 
in the first line of EqJ|. These contributions come from the disclination of the d-vector. 

Using the form of ip(y, z) discussed above, Eq|| I\ and 1% can be obtained as follows. 

l\ = — / dydz 



4 7 [y 2 + (z + R/2) 2 } [y 2 + (z- R/2) 2 } 

R 



7rln 



1 fdydz 1 . 2 ^ 2 ^ 2 /4 



y V ^(y 2 + z2 - R2 /±) 2 + y 2 R 2 J 

= J±)\% (10) 

where ^ is the length scale associated with the spin-orbit coupling defined by £,d(T) = C(T) /firf(T). 

By minimizing f pa i r with respect to R, we obtain the optimal Rq for the lowest free energy configuration of a pair 
of half-quantum vortices and the d-soliton. The optimal Rq is given by 
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Here we have assumed that the £<j > y/eRo/A. Notice that the half-quantum vortices with a d-soliton is possible only 
when K > 1 in order to have Ro > 0. Although we have no information about F", it is most likely that F" ~ 0. The 
ratio between the effective mass and the bare mass, 2*-, is about 4, which means that F\ ~ 9. Therefore, AT > 1 in 
the superconducting state, as one can see from EqJ^. Thus this condition is always satisfied below T c . However, the 
existence of the solution for Rq depends on the value of K. We find that the solution exists only if 1 < K < 1.5. For 
example, for K = 1.5, £,d/R — 0.85. Since the parameter K depends on temperature as shown in Fig. 3, we find that 
a pair of half-quantum vortices with <i-soliton exist only for 0.78 < T/T c < 1. 

Now the free energy of a pair of half-quantum vortices and the d-soltion at the optimal Rq can be obtained as 

'~ - r» c2 ' A1 °? + lJ ^'-2(i^w d + < 12 » 

where A = In ■ In order to examine the stability of the half-quantum vortices, we have to compare f pa ir and the 
free energy of single vortex, /„. The difference is given by 

If f v — fp a i r > for some values of K > 1, a pair of half-quantum vortices are more stable than the conventional 
single vortex. This condition can be rewritten as 

Z\\J 2i K - l )/ 2 {K- l)(*-i)/2 ■ [ > 

Recalling that the solution of Eq. ( ni|) e xists if 1 < K < 1.5, one can investigate the stability condition given by 
One can see from Eq.p] and Eqjl4| that, as long as K > 1, a pair of half-quantum vortices can be stabilized 
over single vortex under certain conditions for the ratio between £<j and A. For example, for K = 1.1 and K = 1.5, 
the conditions for the stability of a pair of half-quantum vortices over single vortex are given by 

^ > 10~ n and ^ > 0.0094 , (15) 
A A 



respectively. One can see that these conditions are easily satisfied. Here we use A/£ = yj (Af,A c ) / (£i,£ c ) = 12.186 which 
is appropriate for Sr 2 Ru04. One can also see that the stability of a pair of half-quantum vortices with the d-soliton 
is determined by the value of K which depends on temperatures as shown in Fig. 3. 

Now let us discuss the relation between I- and d-solitons. It is difficult to estimate the energy of ^-soliton in terms 
of the texture free energy given by EqkL However, it is likely that /-soliton costs much more energy because, if it 
exists, the order parameter given by EqO should vanish inside the Z-soliton. Therefore, if there is a natural passage 
for conversion of Z-solitons to ci-solitons, most of Z-solitons will be converted into d-solitons. 

In summary, assumming that the superconducting state of Sr2Ru04 is characterized by the spin-triplet order 
parameter with broken time reversal symmetry, we investigated the existence of half-quantum vortices and associated 
topological defect; (i-soliton. We showed that a pair of half-quantum vortices attached to a d-soliton can be created 
in the presence of the magnetic field parallel to the a-b plane. It was found that a pair of half-quantum vortices with 
a c?-soliton is more stable than the conventional single vortex for certain temperatures below T c . As in superfluid 
3 He-A, the presence of (i-soliton may be detected as the deficit in the intensity of electron spin resonance signal at 
oj = Qfi. H There should be a clear electron spin resonance signature due to the half-quamtum vortices. Detection of 
the half-quantum vortices by scanning tunneling microscopy (STM) would also provide a convincing evidence for the 
spin-triplet pairing state with time reversal symmetry breaking. 
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auspices of the Department of Energy, supported (in part) by funds provided by the University of California for 
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CAREER award grant No. DMR-9983783 (Y.B.K.) and Alfred P. Sloan Foundation (Y.B.K.). 
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FIG. 1. The spatial configuration of d-vector in the b-c plane given 
by Eq|| The thick line denotes the domain wall with the length of R, 
which is parallel to the c axis. 
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FIG. 2. The spatial configuration of d-vector in the b-c plane given 
by Eq.f?]. The thick line denotes the domain wall with the length of R, 
which is parallel to the b axis. 



7 



K 



1 . 6 



1 . 4 



1 . 2 



1 



0.8 



. 6 




0.7 0.75 0.8 0.85 0.9 0.95 t 



FIG. 3. The parameter K as a function of the reduced temperature 
t = T/T c . 
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